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Abstract 

We study tensor perturbations in a model with inflating Randall-Sundrum-type 
brane embedded in five-dimensional anti-de Sitter (adSs) bulk. In this model, a natural 
in-vacuum of gravitons is the vacuum defined in static adSs frame. We show that this 
vacuum is, in fact, the same as the in- vacuum defined in the frame with de Sitter 
(dS4) slicing, in which the brane is at rest. Thus, during inflation, gravitons on and 
off the brane remain in their vacuum state. We then study graviton creation by the 
brane on which inflation terminates at some moment of time. We mostly consider 
gravitons whose wavelengths at the end of inflation exceed both the horizon size and 
the adSs radius. Creation of these gravitons is dominated by (zero mode)-(zero mode) 
Bogoliubov coefficients and, apart from an overall factor, the spectrum of produced 
gravitons is the same as in four-dimensional theory. "Kaluza-Klein" gravitons are 
also produced, but this effect is subdominant. Graviton spectra at somewhat higher 
momenta are also presented for completeness. 
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1 Introduction and summary 



A special property of brane- world models with infinite extra dimension(s), namely, the 
Randall-Sundrum (RS) model and its generalizations (for a review see, e.g., Ref. [@) is the 
existence of soft bulk gravitons with continuum spectrum of four-dimensional masses. This 
property is of interest in the cosmological context, as in principle it may affect the spectrum 
of long-wavelength perturbations created at inflation, and therefore the large scale structure 
and CMB anisotropy. In inflationary models, perturbations which are of long-wavelengths 
today, are created from zero-point fluctuations of high initial momenta, the corresponding 
"Kaluza-Klein" (KK) zero-point fluctuations are generically not suppressed near the brane 
at early times, and thus may not be negligible. Also, if soft KK gravitons were abundantly 
produced at inflation and survived until recent epoch, they would leave footprints in CMB 
anisotropy. 

The simplest inflationary brane-world model is that of inflating RS brane embedded in 
the five-dimensional anti-de Sitter (adSs) bulk. Naturally, considerable effort has been put 
into the study of the generation of perturbations in this model 0, |], § , in particular, the 
tensor perturbations 0, ||. The latter analyses utilized either brane-based approach || or 
adS/CFT correspondence Q. However, the most straightforward approach to the study 
of the generation of tensor perturbations, which we adopt in this paper, is to solve the 
five-dimensional linearized Einstein equations with appropriate boundary conditions on the 
brane, and calculate the Bogoliubov coefficients. Such a calculation is similar to the study 
of particle creation by moving mirrors: in static adSs coordinates, the three-brane moves || 
and acts as a mirror (with Neumann boundary conditions) for five-dimensional gravitons. 
The calculation of this sort enables one, at least in principle, to evaluate the spectrum of 
created gravitons (not only in the zero-mode final state, as in Ref. |5|], but also in continuum 
states) and also to find the relevant correlators (not only at the de Sitter stage of the brane 
evolution, as in Ref. [[|, but also after inflation). 

In this paper we first consider an inflating brane with de Sitter (dS^ intrinsic geometry, 
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which is a boundary of a part of adS5. We address the following issue. The natural graviton 
vacuum in adSs, cut by the brane, is the state which does not contain five- dimensional 
gravitons moving towards the brane (we refer to this state, somewhat loosely, as adSs- 
vacuum). This definition needs qualification, as it involves the choice of a coordinate frame 
in adSs. The appropriate frame is static, so the brane moves in this frame. From the point of 
view of an observer residing on the brane, another coordinate frame in adSs is appropriate, 
in which the brane is at rest. The latter frame, which corresponds to dS4 slicing of adSs, 
defines its own graviton vacuum, which may or may not coincide with the adSs-vacuum. In 
Section 2 we find that, in fact, the vacuum defined with respect to dS4 slicing is the same as 
adSs-vacuum. This means that the fact that the in-vacuum is defined in static adSs frame 
does not add any physics above and beyond that inherent in dS4 slicing; gravitons on and off 
the brane remain in their vacuum state. This situation is in clear contrast to the case of a 
mirror accelerating in Minkowski space-time. Our result in a sense contradicts expectations 
expressed in Ref. |7j, but we think it is natural. Indeed, consider a brane on which matter 
in its ground state has energy density in excess to the fine-tuned RS value. Then there is 
no reason to suspect that the system would excite itself, i.e., it is natural that quantum 
fields (including gravitons) remain in their ground state, and the expansion of the brane is 
precisely de SitterQ 

Technically, our calculation involves the evaluation of the Bogoliubov coefficients relating 
two sets of graviton wave functions: one set consists of waves moving towards the brane, 
with frequency decomposition defined in the static adSs frame, while another set ||, || |9| 
has fixed frequencies in the frame corresponding to dS4 slicing. We find that the non-trivial 
Bogoliubov coefficients vanish, i.e., the wave functions which are negative frequency in one 
frame, do not have positive-frequency components in another frame. 

In Section 3 we proceed to calculate the graviton creation by the brane on which in- 
flation terminates at some moment of time. We are interested mostly in gravitons whose 
1 In the case of accelerated mirror in Minkowski space-time, there is external force that accelerates the 
mirror. Such external force is absent in our case. 
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wavelengths at the end of inflation are larger than both the CIS4 horizon size H~ l and adSs 
radius We argue in Appendix 1 that for all realistic adSs radii, only this part of the 
spectrum is relevant for the generation of CMB anisotropy at interesting angular scales. 

We find no surprizes. The creation of large wavelength gravitons is dominated by (zero 
mode)-(zero mode) Bogoliubov coefficients and, apart from the overall enhancement factor 
H, the spectrum of produced gravitons is the same as in four-dimensional theories. The 
continuum KK gravitons are also produced, but their number is suppressed. As anticipated 
in Ref. ||, KK part of initial zero-point fluctuations gives negligible contributions into 
production of gravitons in both zero and continuum modes. Hence, the spectrum of large 
wavelength gravitons, created due to inflation on the brane, is indistinguishable from the 
spectrum predicted by four-dimensional theories. 

If the Hubble parameter on the inflating brane, H, is much larger than the inverse 
adSs radius k, there is a range of graviton spatial momenta at the end of inflation, p p h ys , 
in which the gravitons are superhorizon at this moment of time, but their wavelengths are 
smaller than the adS 5 radius, H ^> p p h ys k. For completeness, we consider creation of 
these gravitons in Appendix 4. We find that the spectrum of created zero-mode gravitons 
gets damped at these momenta. One peculiar feature of this case is that the creation of 
KK gravitons in the final state dominates over the production of zero-mode gravitons, in 
contrast to the case of soft gravitons studied in the main text. Still, the energy density 
emitted into the bulk in the part of the spectrum H ^> p p h ys 3> k is a small fraction of the 
energy density on the brane. The emission of gravitons into the bulk is a small effect at 
least for gravitons which are superhorizon by the end of inflation. 



2 CIS4 mirror moving in adSs 
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2.1 in- vacuum in adSs 



The adSs metric in static coordinates is 

ds 2 = -L-(dt 2 - dx 2 - dz 2 ) (1) 
ykzy 

A spatially flat three-brane moving in these coordinates is described by a world surface 

t = t{ V ) , z = z{ V ) (2) 
where rj is the intrinsic time coordinate on the brane. 

Throughout this paper we consider tensor perturbations, i.e., write the perturbed metric 
as follows, 

ds 2 = j]^yl dt2 ~ + h ij )dx i dx j - dz 2 } (3) 

where hij are transverse and traceless. The perturbations are decomposed in spatial Fourier 
modes 

hij{t, x, z) = e 4px e id $(t, z; p) (4) 

where are constant transverse-traceless polarization tensors, and $ obeys the Klein- 
Gordon equation for minimal massless scalar field in adS 5 ^ |5|, 

d 2 3d d 2 



dz^ + -zdz + ^ +P )^ = ° (5) 
With Z2 orbifold symmetry imposed across the brane, the boundary condition for $ is 
Neumann on the brane, 

d n § = (6) 
where d n denotes the derivative normal to the brane. 

The m-vacuum in adSs is defined by choosing the basis of incoming waves (falling 

towards the brane) which have fixed frequencies in the static frame. The solutions to eq. (||), 

which are incoming waves at large z and have fixed frequencies, are0 (up to normalization, 
2 These solutions do not obey the Neumann boundary conditions (|^) on the brane. However, as long as 
one considers asymptotic past, one can use wave packets composed of these functions, and disregard this 
subtlety. 



which is unimportant for our purposes) 



in—adS, m 



e iujt z 2 H^>{mz) 



(7) 



and 



fr H = (d {+) ) 

in—adS, m \rin—adS, mJ 



e- iuJt z 2 Hi 2 \mz) 



(8) 



where Hr, ' are the Hankel functions, the continuous parameter m is positive, m > 0, and 



U) 



\/m 2 + 



(9) 



These solutions behave at large z as follows 



j(±) 

in—adS, m 



3/2 ±i(ut+mz-n/ 4) 



(10) 



As a side remark, we note that the normalization condition is, generally speaking, 



dz y/g g m i{4>* ml d t (t)m - d t (j)* m ,(f)m) = 5 n 



11) 



The factor z 3 ' 2 in eq. © is precisely what is needed to normalize these solutions to delta- 
function 5(m' — m) with the weight y/gg 00 oc z~ 3 . 

2.2 dS 4 slicing 

The frame in which dS4 brane is static is obtained by introducing coordinates rj < and £, 
related to t and z as follows 

t — Tj cosh £ 
2 = — 77 sinh £ 

In these coordinates, the adSs metric is 



ds 2 



{k sinh£) ; 



(d V 2 - rfx 2 ) - de 



(12) 



(13) 



The CIS4 brane is at constant £, 

i.e. in static coordinates it moves along the straight line 

t — Tj cosh £g 

z = — ^sinh^s (14) 

It is clear from eqs. fll2|) and (0) that the induced metric on the brane is indeed de Sitter, 
the Hubble parameter on the brane is 



H = ksinh^B 

and the coordinate rj is the conformal time in dS/t. 

In the new coordinates, the Klein-Gordon equation is 



(15) 



d 2 2d 



d 2 



3 cosh £ d \ , 
+ p 



$ = 



rj drj J rf \ <9£ 2 sinh £ <9£ 
Since the brane is static, the boundary condition (0) takes a simple form, 



(16) 







(17) 



One solution to these equations is the bound state (zero mode), (po(rj), which depends only 
on r). The zero mode which in infinite past has negative frequency with respect to time rj is 



jH = fc f sinhGj— e~ ipr > 



J}- - 

V 



where the normalization factor C\ is determined by eq. (O) written in (rj, £) frame, 

Here k labels the solutions to eqs. (|16"D, (|T7|); for the moment k = k' = 0. One finds from 



eq. pi) 



^(sinhfc) 2 f 

smh £ 



(20) 



The quantity G\ is precisely the same as that introduced in Ref. |J. In the case of slow 
inflation, H <C k (and hence £b *C 1), one has C% = 1, while aX H ^> k one finds Ci = 
vte#/2fc. 

Let us now consider non-zero modes of eg. (|I~6|) . The variables decouple, so the modes 
have the form 

Mv,0 = Mv)-Xk(0 (21) 

where ip K obeys 

* = „ (22) 



drj 2 7] drj rf 

whereas the spatial part of the wave function obeys the boundary condition (|17]) and equa- 
tion 

d 2 3cosh£ d 2 9\ 

+ « + 7 )Xk = (23) 



^<9£ 2 sinh£ <9£ 4, 

Equation ( ^2]) may be reduced to the Bessel equation. Its solution which is negative- 
frequency at large negative r\ is 



^(v) = ^kU^ Hi HfX m ) (24) 



At large negative 77 it behaves as follows 



^ = = jfei ri e -w-i*H . (25) 

V2p 



The normalization factor is chosen in such a way that 

i(Vi +) M^ - c.c.) = 1 (26) 

(where ■?/>i + ' ) = )*)■ The factor k~ 3 r]~ 2 here comes from the measure in (|T§|). 

The spatial wave functions Xk are considered in detail in Appendix 2. Their form is not 
important in this section. It is worth pointing out, however, that there are no bound states 
except for the zero mode, and that the parameter k is continuous and starts from zero. 
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2.3 Vacuum in CIS4 slicing coincides with adSs-vacuum 

Let us see that the wave functions <pk\ which are negative frequency with respect to time rj, 
are also negative frequency with respect to time t. This will mean that the graviton vacua 
defined in the two frames (|1|) and (|12"D, respectively, are in fact the samef]. 

We are going to calculate the Bogoliubov coefficients between the sets (pi ^ and 
^in-adSm = Wn-adS m)* ' These are given by scalar products similar to the left hand sides of 
eqs. ([n]) or (|TJ), but now involving the functions </4~' ) and 4>^_ adSm - We prefer to integrate 
over the hypersurface of fixed rj, and then take the limit 77 — > —00. Hence the Bogoliubov 
coefficients of interest are proportional to 

I sin 



vtt adS oc 



' (sinhOV 



PWVLuJ* - (27) 



where the adSs wave functions are expressed through coordinates rj and £, 

[^-adsj = <Pl ] -a d s, m = e-^ cosh? (^sinh0 2 ^ 2) (-^sinhO (28) 

Let us first note that at large negative rj, both the zero mode 0q ^ and non-zero modes , 
corresponding to dS4 slicing, behave in a similar way, 

<j>P = ve~ ipv X K (0 (29) 



At large negative 77, the integrand in eq. fl27|) rapidly oscillates, and we may use the saddle- 
point approximation. The saddle point £ s is finite, so at large \rj\ we use the asymptotic 
form of the Hankel function iJ^ 2 ^ ( — m? 7 sinh ^) , see also (0). In this way we obtain 

dS-cuiS /T~\ f ^£ /f\ -ioj77 cosh £+im»7 sinh £-ipr? 

v K , m «VM/^ —^ x ^ )e 

\—uj cosh £ + m sinh £ + p + 0(r7 -1 )] (30) 
3 There is a subtlety here. The dS4 slicing (|l^) and correspondingly the wave functions 4> K are defined 
only in the past light cone, z < —t. As long as physics at or near the brane is considered, this subtlety is 
irrelevant. 
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The saddle point is at 



that is 



777 

tanh£ s = — (31) 



Tfl CO 

sinh£ s = — , cosh£ s = — (32) 
p p 

At this point, however, the pre-factor vanishes, up to possible corrections decaying as r/ _1 . 
Thus 

vi S ~ adS = (33) 

as promised. 

A remark is in order. Should we calculate w-coefficients, relating components <f>^ and 
^in-adSm °f the same sig n °f frequencies, the pre-factor would not vanish at the saddle 



point, and we would obtain the result which is finite as rj — ► — oo (the overall factor ■sj \t]\ in 
(p0| ) would cancel H -1 / 2 which would come from the Gaussian integration near the saddle 
point). This may be viewed as a cross check of our calculation. 



3 Gravity waves from inflation on the brane 
3.1 A model for the background 

The main purpose of this paper is to calculate the spectrum of gravitons created by the 
brane which inflates for long period of time, and then enters the Friedmann regime. As long 
as one is interested in gravitons which are superhorizon by the end of inflation, it suffices to 
use the following simplified model for the background. We consider the brane which inflates 
until some moment of time r] < 0, and then does not evolve at all. In a similar model 
in four dimensions, the calculation of graviton creation is very simple, and is presented in 
Appendix 3. Of course, its result coincides with the result of realistic calculations which 
take into account evolution of the Universe after inflation. 
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In static coordinates the trajectory of the brane is the straight line ([14]) for 77 < 770, and 
then another straight line parallel to the t-axis, 




Figure 1: Brane motion in static coordinates: thick solid line represents the world surface of 
the brane, dash-dotted line denotes the past light cone; thin solid (dashed) lines correspond 
to surfaces of fixed r\ (£). 

brane being 

a{z B ) = , 1 1 1 r , (35) 
k\r] \ smh £ B 

Of course, one could rescale rj in such a way that the scale factor on Minkowski brane is 
equal to 1, but it is convenient to keep t]q arbitrary, and make sure, as a cross check, that 
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the results depend on physical quantities. Namely, the physical momentum at the end of 
inflation (and on Minkowski brane) is related to the coordinate momentum p as follows, 

V 

Pphys = — ( — c = pk\r] \ sinh£ B (36) 
a{z B ) 

Comparing this expression with eq. fll5D , we find that the modes which are superhorizon at 
the end of inflation, obey 

(37) 

We will be interested mostly in modes whose physical momentum at the end of inflation 
is small compared to the inverse adSs radius, see Appendix 1. In terms of the coordinate 
momentum, this means 

^i=^ | si nh^«l (38) 
k 

In view of these relations, it will be convenient to consider \t}q\ as a small parameter of the 
problem. 

We note in passing that at H 3> k, the spatial momenta of superhorizon modes may 
be such that the relation ([38]) is violated. Though the corresponding range of momenta, 
H 3> Pphys k, is not particularly interesting in the cosmological context, we consider this 
range in Appendix 4 for completeness. 

3.2 In-out Bogoliubov coefficients: generalities 

We are going to calculate the Bogoliubov coefficients which relate the graviton wave func- 
tions corresponding to in- and out- vacua. The wave functions of the m-vacuum have been 
already described. In our simplified model with Minkowski brane at late times, the wave 
functions of out-vacuum are the standard RS wave functions (modulo factors accounting for 
the fact that the brane is placed at z = Zb)- Namely, the normalized zero mode is 

^ = z B ^ e** (39) 



n 



Non-zero modes are 



where 



^,0 = -S=e ±fcrt ^) ( 4 0) 
y2ui 



< ) - 2 \™L N i( mz B)J2(mz) - Ji(mz B )N 2 {mz) 
~\ 2 y/lN^mzB)]* + [MmzB)]* 
These spatial functions are normalized as follows, 

f°° dz 

j —Pm^m' =S(m- m) (42) 

Jz B Z 

the factor fc 3 / 2 in the expressions (|39|) and ( |40D is due to the measure in (|TT|) . 

Soft KK gravitons have low spatial momentum and frequency (and hence mass m). So, 
we assume in the main text that, together with relations ([J7|) and Q3"5|), their counterparts 

m|r]o| < 1 (43) 

and 

m|?7o| sinh£ B «C 1 (44) 

are satisfied. 

One way to calculate the Bogoliubov coefficients between in- and out- vacua is to specify 
the solution (of fixed frequency) at infinite future, evolve it backwards in time and find its 
decomposition in terms of the wave functions corresponding to m-vacuum. Since negative 
frequency functions in the frame with dS4 slicing are at the same time negative frequency 
functions of adSs vacuum, we are free to use either of these two sets as defining m-vacuum 
(provided we consider the interior of past light cone) 

Thus, the solution of interest coincides in the future (flat brane) region with one of the 
RS wave functions (p. For the solution everywhere in space-time bounded by the brane we 
write 

$ = ip + A$ (45) 

Clearly, the first term here does not contribute to the Bogoliubov w-coefficients, relating in- 
and ottt-functions of different frequencies (the RS wave functions have fixed frequencies in 
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static adS5 coordinates), so the quantity of interest is A$. Using Green's formula, we write 



A<K^£) = -J D adv ( v ,^ v '^ B ) [d?<p(v',Z%'=eB drf 



(46) 



The integration here is performed along the de Sitter brane, i.e. rj' < t]q. The extra term 
A$ concentrates in the past light cone, so we may use the set of the wave functions obtained 
with dS4 slicing, eqs. flTJj) and (|21|), as the set defining the in- vacuum. 

The advanced propagator obeys 



T] 



drf 



2rj-^- + p 2 rf 



where 



d 2 ^cosh^ d 
<9£ 2 sinh £ <9£ 



6 = — - 3 



(47) 



(48 



The form of eq.(f|7|) is chosen in such a way that the Geen's formula holds without any 
additional measure (in particular, the operator <9| enters eq.(|4"T|) without any pre-factor). 

Taking into account that the operator is Hermitean (on wave functions satisfying 
the Neumann condition at £ = with weight sinh~ 3 £ we obtain 



where 



DnM) = ii 0(7/ - vM +) (v)^r\ v ') - ^-\v)^\rf)] i 



k 3 



(49) 



(50) 



Here \n an d are space and time components of the graviton wave functions (including 



the zero mode) in the frame with dS4 slicing, see eqs. and fl2T|) . 
Collecting all factors, we write 

drf 



■ Qirl — rf) x 
sinh 3 £ B r/ 4 V ' /; 



(51) 
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where (f)^ £ ) are the total wave functions in the frame with 0IS4 slicing. These are given 
by expressions (pl|) and (|2T|). 

The advantage of this formula is that the Bogoliubov coefficients follow immediately. For 
example, take <p = (fm ■ Then the f-coefficient, that determines the creation of gravitons 
in this mode from the zero-point fluctuations in the mode <p K , is the coefficient in front of 
0/c + ^(f/, £) in eq. (|5T|) . Namely, we take the limit rj — > —00 and obtain 

*W = "i / Jf^ VM^ fo'tiXv',?)]^* (52) 
or explicitly (and omitting unnecessary primes) 

Similarly, the u-coefficients relating the wave functions with the same signs of frequencies, 
are obtained as follows, 

u m , K = - p ^ff r ^V^fo) I^fo,*)]**, (54) 
/c d smh^B J -00 T 

Thus, calculating the Bogoliubov coefficients amounts to evaluate the wave functions at the 
brane position, Xk{£b), and calculate integrals involving Bessel functions. 



3.3 In-out Bogoliubov coefficients: zero mode to zero mode 

We begin with the calculation of the Bogoliubov coefficients in the case where both initial 
and final wave functions are the zero modes. The final wave function is given by eq. (|31f), 
and in coordinates (??,£) it is 

<pW = -m sinh^e — e ±ipricosh Z (55) 
VP 



whereas the initial wave function is defined in eq. (|Tq) . 

Collecting all factors, we obtain from the representation (EkS 



v 0fi = C lVo r ^(n-Ue-iwa+^s) (56) 



x V 3 \ V 
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For momenta obeying fl37|) and (|38|), the integral in ( [56D is saturated at rj ^ i]q. We find 

^o,o = ^ 57 

which means that at small H/ k, the number of zero-mode gravitons created from fluctuations 
in the zero mode state is 

1 , H 

= 4^f ' ^ = T ^ 1 (58) 

whereas at large H/k one has instead 

3 3 H H 

no '°=8^ SinheB =8^¥' ^B = J» 1 (59) 

In a similar way we get 

"0,0 = -o ( 60 ) 

2pr} 

These expressions coincide with the result of calculation in four dimensional theory (see 



Appendix 3, eqs. (|121|) and ( |122j )), up to overall enhancement factor C\. This is in accord 
with Ref. ||. It is worth noting that the relative phase of m ,o and f ,o coincides too 
(this relative phase enters the correlators of the gravitational perturbations in the post- 
inflationary epoch). 

We note that when calculating the two-point correlators of the gravitational perturba- 
tions after inflation, the overall factor of z^k 3 will appear due to normalization in (|39D . 
This factor may be written as k(Hi] ) 2 = ka 2 . The factor k here accounts for the distinction 
between the five- dimensional field and effective four-dimensional field (i.e., the distinction 
between the five-dimensional gravitational constant and four-dimensional Newton's constant 
whereas the factor a 2 exists also in four- dimensional theory (see Appendix 3, eq.( |117| )). 

We conclude that the (zero mode)-(zero mode) contribution matches the four-dimensio- 
nal calculation, up to the factor C\. 
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3.4 Zero mode to KK modes 



Let us now calculate the Bogoliubov coefficients from initial zero mode to final KK modes, 
which describe the creation of gravitons in continuum from zero-point fluctuations in the 
zero mode of dS4 slicing. Before performing the calculation, we point out that the number 
of created KK gravitons is given by the integral 

dm , 



dm \v m0 \ 



m 



mv m n 



(61) 



Hence, the quantity of interest is 

\/mv mfi (62) 

as it is this quantity that characterizes the number of created particles in a decimal interval 
of KK masses. 



We insert flI5| ) and into expression Q53| ) and obtain 

id 



sinh 2 ^bV^P^ 



'm,0 



(63) 



5 etPV h + ~ ) ( -£? C ° Sh (V sinh £) ] 



(64) 



where 

i\ f d_ 

J\ m \ T V PJ \d£ / 6= t B 
(we changed the integration variable from 77 to —rj). For p and m satisfying fl57D, (P5|), fl4"3] 

and (fHf) , we use the expansion of </9 m at low values of argument, i.e., we insert 



ip m (i] smb. = -z\ 

and obtain, taking into account that zb = \r)a\ sinh^B, that 

uo fm 



(65) 



We get finally 



'm,0 



mv m n 



2p 



sinh d £ B 



-idmy/u . 

sinh £b 



3/2 



(66) 



(67) 



The spectrum rises towards high momenta (recall that t>o,Q oc p 1 corresponds to the "flat 
spectrum"), but the (zero mode)-KK contribution is suppressed, as compared to the (zero 
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mode)-(zero mode) result. In terms of the number of created particles, the suppression 
factor is of order (p\rjo\ sinh^s) 2 or (m|r]o| sinh^g) 2 , which is small in virtue of eqs. fl38|) or 
(|44|). Hence, the creation of KK gravitons from the fluctuations in the zero mode state is 
a subdominant process for soft gravitons. We note in passing that for gravitons of higher 
spatial momenta, the situation is different: we show in Appendix 4 that the creation of KK 
gravitons from the zero mode fluctuations dominates in that case. 



3.5 KK to zero mode 

Let us now consider KK gravitons as initial states and zero mode as the final state. The final 
zero mode is again given by eq. flSoD, while initial KK modes are described in subsection 2.2 
and Appendix 2. 

We are going to calculate the integral (|53|). Using explicit form of the zero mode (|55|). 
we write 



where we again changed the integration variable from rj to —rj. 
Let us begin with the integral over conformal time 



V = k-l^ % ^-\- V ) e*"«*fe (69) 
J\vo\ V 3 

With the explicit form of the time functions given in subsection 2.2, this integral can be 
written as follows, 

W = ^ [ % ^™^e-^H%\ P \ri \ ■ u) (70) 

At small p\f]o\ and p\rjo\ sinh^s this integral is again saturated at its lower limit, u = 1. We 
make use of the representation 

_™ „(i), , | \ e^J iK (p\r] \ -u) - J„ iK (p\r} \ ■ u) 

sinh(7i7cj 



17 



and expand in pr] . We obtain 

7TK/2 



1 



l 0,K 



2 sinh(7rK) (1/2 - m)r(l + m) 
At small k this expression behaves as 



P\Vo\ 



+ (k 



(72) 



2i 1 



<0,K 



sin 



7T ft 



ft log 



P\Vo\ 



K < 1 



(73) 



while at large k it behaves as 



T _ 1 n mlQg(p|T ?0 |/2)+i«(l- 

- 1 0,K 



log/c) 



, « > 1 



(74) 



It remains to use the formulas for spatial functions X«(£b) to obtain the Bogoliubov 
coefficients. Let us begin with the case sinh^s = H/k <C 1. The contributions to the 
number of created zero-mode gravitons, coming from zero-point fluctuations with small and 
large k, have different form. At k <C sinh -1 we use eq. (|107|) and obtain 

1 



v ,k = ~ 



V2 



- + i- 



4^ L 2 



^p\t]q\ sinh VretaiiriK • Io, K , ft^Csinh -1 ^ (75) 
At k ^> sinh" 1 £b, equation (|113|) with cos# K = 1 is appropriate, and we find 

-t=Vp\Vo\ sinh^ B • h, K (76) 



Making use of eqs. ([73] ) and ([71"D we find that the integral for the number of created gravitons 
is saturated at k <C sinh -1 and is of order 



dn \vq, k \ 2 ~ p|^o| sinh 2 ^ , sinh^ B <ti 1 



(77) 



Clearly, it is suppressed as compared to the (zero mode)-(zero mode) contribution to the 
number of created gravitons, eq. (]58|) , by a factor (p|?7o|) 3 sinh 2 ^ = (p p h y s/k) 2 ■ (p p h ys /H). 



Let us now turn to the case sinh£ B = H/k ^> 1. In this case we use eq. flll3|) and find 



v 



,k ~ y/p\Vo\ sinh^ B • J 0(t cos 6 K 



(7E 
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For k < 1 and k 1, one has cos8 K = 2k/3 and cos# K = 1, respectively. Making use of 
eqs. ([T3|) and ( |74"|) we find that the integral for the number of created gravitons is of order 

J dn \v 0jK \ 2 ~ p\r] \ sinh^ B , sinh^e > 1 (79) 

Again, it is suppressed as compared to the (zero mode)-(zero mode) contribution, eq. (|59|) . 
but now by a factor (p|?7o|) 3 = (p P h ys /H) 3 . 

We conclude that the contribution of the zero-point fluctuations of the KK modes into 
the production of soft gravitons in the zero mode state is negligible. 

3.6 Continuum to continuum 

Finally, we calculate the Bogoliubov coefficients in the case where both initial and final wave 
functions belong to coninuous spectra. We are again interested in the quantity ^pmv m ^ K (see 
the discussion in the beginning of subsection 3.4). 

The expression (|5^) now reads 

*W = 4= r % e-™ /2H %(m) ^[e— h W^inhOk=, s (80) 

We again use the expression (^) and obtain 

V^v m>K = -- ■ \vo\ 3/2 X k (Zb) m^J J™ ^ e Mv \u C os^ Be -^/2 H (i) {plml . u) (gl) 

The integral here is the same as in eq. (|70|) , and it is dominated by the lower limit of the 
integration. Comparing eq. fl8T|) with eqs. (|68|) and ([7(]) we find that at small p and m, the 
creation of KK gravitons from KK zero-point fluctuations is further suppressed as compared 
to the case of subsection 3.5. For the number of created gravitons, the suppression factor 
is of order (p\r] \ sinhy^) 2 or (m|?7o| smliy^) 2 . Thus, the "KK to KK" contribution to the 
creation of soft gravitons is completely negligible. 
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3.7 Discussion 



We have found that the main source of created gravitons is the initial zero-point fluctuations 
in the zero-mode state existing in the frame with dS4 slicing. The contributions of initial 
KK modes are negligible in all cases. As shown in Appendix 4, this general property holds 
not only for soft gravitons considered in this section, but also for gravitons of somewhat 
higher spatial momenta. 



Soft gravitons obeying ( p7|) and (|38|) (and also (|4^) and fl44| ) for KK modes) are created 
predominantly in the final zero-mode state. This is not, however, the general property of this 
system: we find in Appendix 4 that the situation is inverse for gravitons whose wavelengths 
at the end of inflation are shorter than the adSs radius. For the moment, the novel feature of 
the inflationary brane-world model — the creation of KK gravitons, or, in other words, the 
emission of five-dimensional gravitons into the bulk — appears of academic interest only. 



The authors are indebted to N. Arkani-Hamed, P. Binetryi, N. Deruelle, S. Dubovsky, 
D. Langlois and P. Tinyakov for helpful discussions. This work has been supported in part 
by RFBR grant 99-02-18410, CPG and SSLSS grant 00-15-96626, CRDF grant (award RP1- 
2103), Swiss Science Foundation grant 7SUPJ062239. The work of D.G. was supported in 
part also under RFBR grant 01-02-06035. The work of S.S. was supported in part also 
under RFBR grant 01-02-06034. 

Appendix 1. Smallness of graviton momenta. 

We consider in the main text the creation of gravitons whose physical momentum at 
the end of inflation is smaller than the inverse adSs radius, 

Pphys < k (82) 

Here we argue that only these gravitons are of interest from the point of view of CMB 
anisotropy. For this purpose, we estimate the minimum possible present momentum p of 
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gravitons which at the end of inflation have 



Pphys k 

We will see that in any feasible scenario, the present wavelength 1/po is too short to be 
relevant for CMB anisotropy. 

In the first place, we are dealing with gravitons which are superhorizon by the end of 
inflation, 

Pphys <C H 

Thus, we only have to consider the case 

#» A; 

There is nothing wrong with this case, but the cosmological evolution is somewhat different 
from the conventional four-dimensional cosmology. As shown in Ref. |TIJ, instead of the 
conventional Friedmann equation one has in this case 

H 2 = ■ 4P 2 (83) 

where n\ = M^Jj is the five-dimensional gravitational constant. With = kMp t this 
means 

PeA . ~ HkM 2 Pl (84) 

This determines the energy density just after inflation (hence the subscript, meaning "end 
of inflation"). 

Let us set an upper limit on the redshift from the end of inflation to today. Consider 
an extreme case of very long post-inflationary preheating, at which the equation of state 
corresponds to pressureless matter. Let us assume, again as an extreme case, that this 
regime terminates only at the nucleosynthesis epoch, Tjss ~ 10 MeV. Then 



rp / \ 1/3 

do Clo dNS J-NS I Pe.i. 



dNS CLe.i. To V^TVS 
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15) 



where the subscript refers to today. 

Now, the Hubble parameter at inflation, H, must be smaller than the fundamental scale 
M (5) , i.e., 

H < (kM 2 Pl ) 1/3 (86) 

In fact, a stronger bound comes from gravity waves themselves, as their amplitude (coming 
from (zero mode) - (zero mode) Bogoliubov coefficients, with account of C\ ~ y^H/k) is 

l ~ikJi (87) 

It must be smaller than about 10~ 3 , so instead of (B6|) we have 

H < (5 2 kM 2 Pi y/ 3 (88) 

where 5 ~ 10 -3 . 

For the present wavelength of a mode whose wave number just after inflation was equal 
to k we have 

p k k T V t ns J 
The smaller k, the larger l , even if one takes into account (|88|). For k' 1 ~ 1 mm one obtains 
from (§§ 

H < 3 ■ 10 6 GeV (90) 

and then from (^) one finds 

l < 10 23 cm ~ 30 kpc (91) 

As claimed, this is too short scale to be relevant for the CMB anisotropy. 

In all above estimates we pushed everything to extreme, in more realistic cases the scale 
/o should be much smaller. 

Appendix 2. Graviton modes in the frame with de 
Sitter slicing 
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Small sinh^s and ftsinh£g 

In terms of the variable u = cosh 2 £, eq.(^) is the hypergeometric equation 

u(l - u)d 2 uX + Q + d uX - + ^ X = (92) 

(the subscript in Xk is omitted). Another way to cast eq.(|23|) into hypergeometric form is 
to introduce the variable v = — sinh 2 £. Then one has 

v(l - v)dlx + (-1 + \v\ d vX - (1 + ^pj X = (93) 

It is then convenient to choose two linear independent solutions as follows. The first solution 
is 

X i = F(a,/3;i;cosh 2 ^ (94) 

where 

3 , Ac ^ 3./^ /rt _ N 
a = /3 = i- (95) 

4 2 ' P 4 2 v ; 

and the hypergeometric function F should be evaluated at a fixed side of the branch cutQ 
extending from 1 to oo (say, above the real axis). This form of the solution is clear from the 
representation (|9~2"1). The form of the second solution comes from the representation (|93|) 

X 2 = sinh 4 e F Q - or, - - /?; 3; - sinh 2 £j (96) 

The reason to choose this particular pair of solutions is that they behave quite differently 
at small £: the first one behaves as 

T(l/2) 

Xl = r(l/2-a)r(l/2-/3) (1 + a ^ sinh ^ + 0(sinh 4 e log(sinhO)) (97) 
whereas the second one behaves as 

X2 = sinh 4 e (98) 
As a by-product, these formulae tell that xi an d Xi are indeed linear independent. 



4 Another choice of the side of the branch cut would correspond to a linear combination of the original 
solution and the second one, X2- 
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We now write the required solution in the form 



X = A lX i + A 2X 2 (99) 

and impose the Neumann condition fll7D . We first perform the calculations assuming that 

sinh£B<Cl and Ksinh£ B < 1 (100) 

The opposite case is easy and will be considered later. In the case (|100| ) we use the small-£ 
expansions (|97|) and ( p8|) near t; = (,b and obtain 



l2 2r(i/2 - a)r(i/2 - 3) , v 

Note that A\ is small compared to A 2 for small sinh£ B . 

We now wish to calculate the normalization factor. The following normalization of 
continuum modes is appropriate, 

~^rXx«' = S ( K - K ') ( 102 ) 

As is usual for continuum spectrum, the normalization factor is determined by large-£ 
asymptotics of the solution. It is straightforward to find the asymptotics of the first solution, 



Yl - r(i/2)r(-iK) (coah a*-* e * 

Xl " r(-3/4-m/2)r(5/4-m/2) (COSh ^ 6 

+ r(i/2)r(zK) {oodk rf«* c * 

+ r(-3/4 + m/2)r(5/4 + m/2) 1 U 
where phases 6 1 and 6' are not important for our purposes. The asymptotics of the second 
solution is 



r(3)r(-m) . 

X2 ~ r(5/4 - m/2)r(7/4 - w/2) ^ 



IY3)r(/,v) (sinhO^ +iK (103) 



3 



T(5/4 + m/2)r(7/4 + m/2) 
At large £ one neglects A\Xi term in (|99|) . Making use of (|103|) , one finds 

r(5/4 + m/2)r(7/4 + m/2) 



.4, 



2^ 



r(3)r(«) 
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(104) 



(sign here is chosen for future convenience). This expression can be simplified by making 
use of the identities, 



. , 5 k\ _ (1 K 
' '3 + *2 r 4 + *2 



1 k\ ( 3 K . , 

4 +i 2 r 1 -4 +i 2 1 ' 11 



3 K 
-7+1- 

4 2 



2tt 5 



cosh (jr k) 



\TUk) 



1 K 
- +1- 

4 2 

7T 



3 k 
- + i- 

4 2 



k sinh (7m) 



We find 



1 /k , 
■~\ — tanhK 
2V 2 



1 K 

4 +! 2 



3 /c 
4 + *2 



(105) 



(106) 



For the calculation of the Bogoliubov coefficients we need the value of the wave function at 
the brane position, £ = £ B . Since A 2 X2(£,b) ~ v4 2 sinh 4 £ B and AiXi(£b) ~ A 2 k~ 2 sinh 2 
the part is dominant on the brane. Making use of eqs. (|97D and ( |101| ) we find 



B 



-A, 



I- 4- 7-I 2 

4 ^ *2 I 



sinh 2 £i 



1 II + i*\ 

k tanh k\4 • sinh 2 £ 



B 



(107) 



We recall that these formulae are valid provided the relations ( |100| ) are satisfied. We now 
turn to the opposite case. 



Large sinh^B and/or ft sinh 

In the opposite case, 

sinh<^B 3> 1 and/or /tsinh^ 1 
it is convenient to introduce the wave function 

X(0=X(0 sintT 3 / 2 f 



and cast equation fl23|) in the form of the Schrodinger equation 

15 



-d + 



4 sinh £ 
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(108) 



(109) 



(110) 



In the case of interest, the potential term in eg. ( |110|) may be neglected, and the solution 
(normalized to delta-function) becomes 



-^sinh 3/2 £ cos[k(£ - £b) + 0*] (Hi; 



X 

V 71 " 

The phase here is determined from the Neumann boundary condition on the brane, eq.fll 
One finds 

3 cosh^ B 

tan ^ = ^ — r ~TT~ 112 

2 k smh £b 

Thus, the value of the wave function at the brane position is 

Xk {£b) = ^sinh 3/2 ^cos0 K (113) 



7T 



Note that for parameters obeying eq. ( 108 ), the value of the wave function on the brane is 
always smaller than sinh 3 ^ 2 ^- 

Appendix 3. Bogoliubov coefficients in four-dimen- 
sional theory. 

In the four-dimensional theory one considers the metric which is de Sitter at inflationary 

stage 



and flat afterwards, 



ds 2 = -jj^jy( d V 2 ~ dx 2 ) , v<Vo (H4) 



ds 2 = {drf - dx 2 ) , t] > Vo (115) 



? 

Properly normalized solutions in the de Sitter and flat regions are 



V%P V P 

and 



±) = -^Le ±4p " ( rj ± - ) (116) 



= -iL^e"™ (117) 
V2p 



respectively 
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In the de Sitter region, one writes the solution that becomes <// ) at late times as a 
linear combination 

$ = u * ft-) + v< j)(+) (118) 

and finds the Bogoliubov coefficients by matching this solution and its derivative to at 
r] = r] . In this way one obtains 

u=e -2ipuo_L_ ( U9 ) 



2m 





u = 1 — (120) 

At small p?7o (physical momentum just after inflation small comared to H), these become 

v = ^~ (121) 
u = — (122) 



This is the standard "flat spectrum" . 

Appendix 4. Creation of gravitons of higher mo- 
menta 

At H >> k, i.e. simile 1, there exists a range of momenta in which the gravitons are 
superhorizon at the end of inflation, 

p\Vo\ < 1 

and yet their physical momentum at that time exceeds the inverse adSs radius, 

p\r] \ sinh^B > 1 (123) 

Here we present for completeness the results of calculations of the Bogoliubov coefficients 
in this case. For the masses of final KK gravitons we assume also 

m\i] \ < 1 , m|?7o| simile > 1 (124) 

so that their physical frequency is small compared to H and large compared to k at the end 
of inflation. 
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Zero mode to zero mode 

The integral flSBQ is equal to 



ft pimp cosh § B 

%o = — (125) 

pi] prjocosh^B 

where, as before, C\ = a/3 sinh £b/2 at large sinh^- Comparing this result with eq. ( |57j ) 
and recalling (|123|) , we see that the spectrum of created gravitons is damped at higher 
momenta due to the extra factor (pr] Q cosh£ B )~ 2 appearing in ra = |t>o,o| 2 - 

Zero mode to KK modes 



To evaluate the integral (|6^), we use the approximation 



P) , ,„ p«*«7 coshes 



[e^ cosh V m (r/sinhO] = (?? sinh£ B ) 5/2 



<B V 71 " 



x (zcj cos (m(|?7o| — v) sm h£B) + m coth£ B sin (m(|?7o| —rf) sinh£ B )) (126) 

We obtain 

™m,o = -^3= cotHse^^ - Jl (127) 



2v / 2~7T ?7o(cu 2 coth m 

This result is to be compared with eq. fl5"7p. We see that the spectrum no longer increases 
towards high momenta; instead, it turns down at large p and m. It is of interest also to 
compare eqs. ( 127|) and ( |125| ) and find that for higher momentum gravitons 

Kol 2 i 



\Vmv mfl \ 2 p\t] \ cosh^ B 
This means that as far as the zero-mode initial fluctuations are concerned, the production 
of KK gravitons dominates over the production of the zero-mode gravitons in the range of 
momenta considered in this Appendix. This is in contrast to the case studied in the main 
text. 

KK to zero mode 

The Bogoliubov coefficient reads 

smh t, B 
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where is given by Eq. (|113|) and 



v 7F 

'0,K 



l°° % ^ 0lco8HB - u e- J fH^(p\r l0 \-u) 



2 

This integral is again saturated at the lower limit of integration, so in the integrand we 
substitute the expression for the Hankel function at small value of its argument. In this way 
we obtain 

Jtx e ipM cosh fs 1 / //>|//o| V" * ,_J " 

Iq.k = i—: r-n — \ — ; — ; r^— — ^— / 1 i ~ \ - (« -> -«) 



2 sinh(7rK) p|r/ |cosh^ B I \ 2 / r(l + i«) 
The final estimate for the number of created gravitons is 

J dn \v , K \ 2 ~ (p\rj \ sinh^) -1 (128) 

Comparing this estimate with eq. ( [79] ) we find that the spectrum is again damped at higher 
momenta: the ratio of eqs. (|128|) and (]79f) is of order (p\i] \ sinh^)™ 2 . Also, it follows from 
eqs. ( |125| ) and ( |128] ) that the production of the zero mode gravitons from KK modes is 
suppressed by the factor (p|??o|) 3 m comparison with the production of zero mode gravitons 
from the fluctuations in the zero mode. This suppression is the same as for soft gravitons 
considered in subsection 3.5. 

Continuum to continuum 



The Bogoliubov coefficients are given by eq. (|80|) where Xk(£b) is determined by 



eq. (|113|) , and the approximation ( |126| ) can be used. Thus, one has to calculate the in- 



tegral 

It? I p~ 7rK / 2 r°° 

v m , K = — sinh£ B cos9 K —=r- / du H^(p\r) Q \ ■ u ) e Hwl«*he fl .ti 

x (iuj cos (m|r/o|(l — u) sinh^s) + m coth£s sin (m|r/o|(l — u) sinh^)) 

It is saturated again at the lower limit of integration. We find 

1 cos# K e iLj|rw|coshfs 



x 



2 \>2uj sinh7m 

uj — m coth£# 
uj coth£ B — m 



e W2 f P \ Vo \\ m f 1 _ K 



r(l + i«) \ 2 J \ \rj \(Ltj coshes — msmh.£ B ] 



— ( K — > — K 



+ (m -> -m) } (129) 
29 



One can show that this expression leads to the number of created gravitons suppressed by a 
factor {p\T]o\ sinh^s) -3 as compared to the case of soft gravitons, see eq. flSip. Furthermore, 
comparing eqs. ( |129| ) and ( |127| ), one finds that the production of KK modes of higher 
spatial momenta from KK modes is suppressed by a factor Qo|?7o|) 3 in comparison with the 
production of the KK modes from the zero mode. 

Discussion 

For momenta obeying (|123|), the dominant effect is the creation of KK gravitons from 



initial fluctuations in the zero mode. In five-dimensional language, this corresponds to the 
emission of the five-dimensional gravitons into the bulk. Let us estimate the energy per unit 
three-volume, emitted into the bulk in the form of superhorizon gravitons, 

1 f dm d 3 p — 2 

(^superhorizon ~~a"/ \ / 77i vT ^ V ^^m,0 

where the factor a -4 (zb) accounts for the fact that p and u are conformal, rather than 
physical, three-momentum and frequency. Making use of eq. (|127|) , and setting the upper 



limit of integration equal to p ma x ~ rrimax ~ |?7o| 1 , we obtain an estimate 



1 

a 4 {z B )r]$ 



^superhorizon ~ „AI ~ \ 4 ~ ^ 



In view of eqs. (|8^ ) and (|88D , we find 



^superhorizon ^ ^2 
Pe.i. 

This means that only a small fraction of energy leaks from the brane into the bulk in the 
form of gravitons with momenta smaller than H. The emission into the bulk is a small 
effect, at least for superhorizon modes. 
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